Stable Magnetic Universes Revisited 



O 

(N 



5-H 



T. TahamtaiE and M. Halilsojlll 
Department of Physics, Eastern Mediterranean University, 
G. Magusa, North Cyprus, Mersm 10 - Turkey. 

A regular class of static, cylindrically symmetric pure magnetic field metrics is rederived in a 
different metric ansatz in all dimensions. Radial, time dependent perturbations show that for 
dimensions d > 3 such spacetimes are stable at both near r ~ and large radius r — >■ oo. In a different 
gauge these stability analysis and similar results were known beforehand. For d — 3, however, 
simultaneous stability requirement at both, near and far radial distances can not be reconciled 
for time - dependent perturbations. Restricted, numerical geodesies for neutral particles reveal a 
confinement around the center in the polar plane. Charged, time-like geodesies for d = 4 on the 
other hand are shown numerically to run toward infinity. 



I. INTRODUCTION 



In 4— dimensional spacetimes the electric - magnetic duality symmetry of the Maxwell equations is an important 
property which can not be satisfied in other dimensions unless different form fields other than 2— forms are introduced. 
\ For this reason a dyonic solution admits a meaningful interpretation only in c? = 4. The Reissner - Nordestrom (RN) 
solution constitutes in this regard the best example which has both electric and magnetic solutions in a symmetric 
manner. In other dimensions ( d ^ 4 ) similar duality properties can, in principle, be defined as well but physical 
interpretation corresponding to electric and magnetic fields turn out to be rather abstract. For such reasons, in order 



I ' to avoid complications due to the absence of a tangible duality, pure electric or pure magnetic solutions seemed to 



attract considerable attention. This amounts to only half of the Maxwell equations, the other half being trivially 
satisfied. From this token we wish to resort here to the pure magnetic solutions which yield a completely solvable 
class without much effort. From the physical side, occurrence of pure and very strong magnetic fields associated with 

■ astronomical objects such as magnetars motivate us to search for such solutions in general relativity. 

^ ' It was Melvin, who first studied such cylindrically symmetric parallel magnetic lines of force remaining in equilibrium 
under their mutual gravitational attraction in d = 4[1]. Later on, generalized version of the Melvin's magnetic universe 
' was also considered [2[. The Melvin universe is invariant under rotation and translation along the Zi axis orthogonal to 
: the polar plane (r, ip). Thorne popularized the Melvin universe further by showing its absolute stability against small 
' radial perturbations [Sj]. Due to this stability property it can be presumed that astrophysical objects emitting strong 
I beams of magnetic fields may everlast in an accelerating universe. Additionally, in d = 3 [J and 5 dimensional 0| 

■ cases also pure magnetic field solutions were found and their energy content investigated Q . Pure magnetic solutions 
T-H ^ in higher dimensions are also known to exist in string. Lovelock, Yang- Mills, Born-Infeld and other theories 01 • 

^ . In this paper, we present in a particular cylindrically symmetric metric ansatz, a class of non-singular, source-free, 
• ^ ' static, pure magnetic solutions to Einstein-Maxwell (EM) equations in all dimensions. In a different metric ansatz 
^ , these solutions were known previously Our principal aim is to investigate the stability of such magnetic universes 

^ ' against time dependent small radial perturbations and explore the possible role of dimensionality of spacetime in such 
- - - matters. It has been known for a long time that for d > 3 these kind of magnetic solutions are all stable 8] . We verify 
these results once more in a different metric (i.e. non canonical) ansatz with supplement of the d = 3 case.. We show 
that for d > 3 the metrics are stable against small perturbations at both near axis Zi and at far distance away from 
Zi. We observe also that when d = 3 these two regions behave differently. Namely, the metric can be made stable at 
r « or, at r — ^ oo, but not simultaneously, which we phrase as 'weakly' stable. The solutions justify once more the 
impossibility of cylindrical magnetic field lines implosion and therefore formation of such black holes. We investigate 
the time-like ( d > 4, with fixed polar angle) and null ( d = 3,4 ) geodesies for neutral particles numerically. Only for 
d = 4 and 5 we were able to obtain exact integrals, albeit in non-invertible forms, of the geodesies equation. In each 
case a confinement of geodesies is observed to take place near the central region. Due to its physical importance we 
consider also the time-like geodesies of a charged particle. It turns out that such geodesies can not be confined and 
in their proper time they diverge to infinity. 
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Organization of the paper is as follows. In Section II we present our metric, field equations and solve them in 
d— dimensions. Perturbation analysis of our system follows in Section III. Geodesies motion is studied in Section IV. 
Our results are summarized in Conclusion which appears in Section V. 

II. METRIC AND SOLUTIONS OF FIELD EQUATIONS IN d-DIMENSIONS 

Our d— dimensional static, cylindrically symmetric line element ansatz is given by 

ds^ = f (r) (^dt' - dr^ - Edzf^ " 7^^^'' d>3 (1) 

in which / (r) is a function of r to be found and bo and k are constant parameters. Also the pure magnetic 2— form 
field is chosen to be 

F =Fr^dr A dip, (2) 
where Fr,p is the only non-zero component of the electromagnetic field. The energy momentum tensor is defined by 



Ai^Ti = -F.kF^'' + J(5jF„„F™" (3) 



which admits the non-zero components 



Ti = diag [TO = T| = T| = . . . = -T^ = -Ttl] = ^^r^F''''- (4) 

We note that our choice of indices {0, 1,2, . . . ,{d — 1)} denote {t,r,zi, - ■ ■ , Zd-z, (f} and the energy conditions satisfied 
by this energy-momentum tensor are discussed in the Appendix A. Prom the Einstein equations, T- = G\ Eq. (4) 
implies that 

<Jo — — *j^3 — — — -^d-v y^i 

Prom these relations we can write 



where 



G\ 



G" 



d-l 



(6) 
(7) 
(8) 



Gl 



G\ = 



1 



Gt 



Apr 
1 

Jpr _ 
1 



2// "r (fc - (d - 3)) + rf ^ i^-k^ + {d - 6) k - {d 
rf^({d-2)k-{d- 2) fc^") - 2ff{d - 2) 



Apr 



-2ff "r{d-2)-rf'^{d-2) 



id -7) 



Prom equation (6) we obtain the differential equation 

2rf "/ + (fc_2)rf 2-2//'=0. 

whose solution for d > 4 is 

/(r) = {kr^+Ci)i. 



(d-3) 



2ff{2k-{d-4)) 



(9) 
(10) 

(11) 

(12) 

(13) 
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By putting this result into Eq. (7) for finding k we find out that for d > 4, fc is d — 3. On the other hand, for d = 3 
the solution turns out to be 

/(r) = C2e^^" (14) 

for the integration constants C2 and c. For convenience we make the choices Ci — d — 3 and C2 = 1, so that the 
solution can be expressed by 

f(d-3)^(r; + l)A d>4 
Accordingly, our line element takes the form 

rf,2 ^ I {d-3)^{r^ + 1)^ (dt^ dr- gdz,f) - (,_3-;S^,). V, d>4 ^^^^ 

\ e^"-" {de - dr^) -rhld^^, d = 3 

We note that these solutions are not new, for they coincide with those of @ (for d > 4) and for (d = 3), respectively. 
It can easily be seen that for r — > it reduces to the following form 

(d-3)^{de-dr^- Zdzf)-j0^dip^ d>4 ^^^^ 

{dt^ - dr^) -~r%ld^^, d = 3 

52 

This represents a conical geometry signalling the existence of a cosmic string near r = 0. By choosing (^ii^^yi = 1 for 
d > 4 and 60 = 1 for d = 3 , we have Minkowskian metrics as one approaches the axes Zi. The solution (16) is a 
singularity free magnetic universe in d— dimensions in analogy with the Melvin space time. 
From Maxwell's Eq. it fallows that 



So 



d > 4 



I d = 3 



{Bq = an integration constant) 
which implies that the magnetic field behaviors as a function of r are 



^'■'^ 1 r d-3 



The marked distinction between d = 3 and d > 4 cases can already be seen from these behaviors. Accordingly the 
energy density reads 



2d-4 2d-4 



d> 4 



Tq" = <( 87r(d-3)^J^(l+r2)T^ (19) 



^e-""" d = 3 



while the Ricci scalar for the metric (16) is 



l-d 

-4(d-4)(d-3)^ 



d> 4 



2ce-=''' d = 3 



Similarly, the Kretchmann scalar has the behavior 



d > 4 



4(d-2) 4(d-2) 

K - <^ (rf-3) <i-3 . (21) 

[ 12c2e-2'='- d = 3 

It is observed that regularity at r — >■ oo dictates us to make the choice c > for the integration constant. 
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III. PERTURBATION ANALYSIS 



In this section we perturb the metric and magnetic potentiaL Since the case d = 3 forms a special case we consider 
it separately. Similar analysis was carried out by Gibbons and Wiltshire Q where they used the canonical metric 
ansatz. We shall show below that their results can also be obtained in a different metric ansatz. 



Our line element is 



1. The case for d — 3 



ds^^f{r,t) (df-dr^) 



9{r,t) 



where 



and the magnetic potential is expressed by 



The unperturbed functions are 



f{r,t) - fo{r)+eu{r,t) 
g{r,t) = go{r) + ew{r,t) 



A^{r,t) ^ A^{r) +ea{r,t). 



(22) 



(23) 



(24) 



fo{r) - e= 
9o{r) = 1 



Bobor^ 



(25) 
(26) 

(27) 



where u (r, t) , w (r, t) and a (r, t) are the perturbed functions. Since e is a small parameter we assume that w in 
our analysis. We use Einstein's equations to find the perturbed functions. The differential equations satisfied by the 
perturbed functions are 



dr 



d^u(r,t) d'^u(r,t) , , du(r,t) . -i ■) , n du(r,t) 



dr 



dr 



(28) 



d'^wir.t) dw(r,t) d^w(r,t) 
1 TT^ K 2 — r -—: — — 



dr^ 



dr 



d'^a{r,t) d'^a{r,t) lda{r,t) B^hor dw{r,t) 



Q^2 J, Qj. 

This system of differential equations admits the solutions 



dr 



= 0. 



(29) 
(30) 



u{r,t) 



2 f ArBn 

e ( [Bih{ar) + B2Ki{ar)] + 



V ^0 



(31) 



ere 



( — ar — a) 



{-(sinhar + coshQ;r)(£;i + E2)e^°"" + e^"(sinhQ;r - coshar)(-i;i + £'2)} 



{{crE2 — aEi) coshar + {crE\ — aE2) sinhar}) 



(32) 
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a{r,t) = re-"* (Bihiar) + B2Ki{ar)- (33) 
^j^e-"''-°' {{sinhar + coshar){Ei + £■2)6^°''' + e'^^'isinhar - coshar){-Ei + £2)}^ , 

where Bi, B2, Ei, E2 and a are all integration constants while Ii{ar) and Ki{ar) are the modified Bessel functions 
of order one. At all times t, to a first order in e , the locally flat nature of the metric near r = will not be altered 
This implies when r — >■ 



(—r^^^Ki {ar) - crEie-""-"' {(sinhar + coshar)?.^'^" + e2"(sinh ar - coshar)} (34) 



+Ei {cr sinh ar — a cosh ar}) , 



, -atj^ sinh ar 

w{r,t) — e El , (35) 

r 

a{r,t) =re~"* (^B2Ki{ar) - ^^^£;ie""''"" {(sinha + cosha)e2"'~ - e2"(sinha - cosha)}^ . (36) 

It is seen that we can choose Bi = and i?2 = to have finite limit when r — >■ 0. However, for r — > cx) it can 
be checked that w{r,t) diverges and since go{r) — 1 this implies that the ratio of ^^^jyy grows indefinitely. No other 
choice of constants suffice to eliminate this divergence. As a result the perturbation converges for r — > but diverges 
for r — )■ cx) in d = 3 case. Let us note that for the time - independent perturbation both, near and far - region 
perturbation terms become convergent. 

2. The case for d = 5 and d> 5 
Our perturbed line element now for d = 5 is of the form 

J.2U2 

ds^ = f(r, t) (dt^ - dr^ - dx"^ - dy^) - —-A-difi^ (37) 

g{r,t) 

From Eq. (23) and (24) we introduce in analogy, the perturbed functions in which the unperturbed functions are 
given by 

/o(r) = 2(r2 + l) (38) 
.go(r) = 4(r2 + l)2 (39) 

Bobo l_ 

■ (r2 + 1) 

From the Einstein's equations we obtain to the first order in e the following differential equations for u{r, t) and w{r, t) 



Kir) = (40) 



df^ 2 at2 

Integration of these equations yield 



.2!^_2i^_„„,t) = o (41) 

2(^" + l)^-i^^0 (42, 



u{r,t) = e-"^* [r{ph{ar) + qKi{ar))] (43) 
w{r,t) = Ae-'^*[r{ph{ar) + qKi{ar))]{r^ + 1) (44) 
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for integration constants p^q and a, and Ii[ar) and Ki{ar) are the modified Bessel functions of order one. At all 
times t. to a first order in e , the locally flat nature of the metric near r = remains intact. Also, for all times to 
a first order in e , the static metric for r — )• oo will not be altered. In this case we see that perturbed functions are 
flnite at r = 0, and when r ^ oo they go to zero. 

A similar analysis has been carried out for cZ = 6, 7, ... so that we found the general solution for the metric functions. 
In each case we have obtained the following relations between w (r, t) and u (r, t) 



w{r,t) = {d-Z) 

/o(r) = (d-3)A(r2 + l)^ 
g,{r) = {d-^)\r-' + lf. 

Solutions for w (r, t) and u (r, t) are given in all dimensions as follow 



(45) 

(46) 
(47) 



w{r,t) = (d-3) ''-^ {r^ + l)^r=3u{r,t) , 
u{r,t) = e-"^d - 3)^ {ph {ar) + qKi (ar)) 



(r2 + l)< 



(48) 
(49) 



in which p and q arc constants. These are both finite at r — > and the ratio of them when r — > oo go to zero. 
Now, from Maxwell's equations we attempt to find the perturbative solution for the magnetic potential. We take the 
potential in the form 



A^{r,t) = A^{r) + ea{r,t) 



where 



2(d-3)2 (r2 + l)' 



The difi'erential equation satisfied by a (r, t) becomes 



d^a (r, t) _^ d^a (r, t) _^ (Sr^ - l) da (r, t) 



boBor 



9*2 



dr^ r (r^ + 1) dr 



u (r, t) 
(r2 + 1)^ 



in which 



u (r, t) = e""' {d - 3) ''-^ (ph (ar) + qKi (ar)) 



(r-2 + 1)3=3 



(50) 



(51) 



(52) 



Upon substitution for u (r, t) we see that for all dimensions the equation satisfied by a (r, t) takes the form 



d'^a (r, t) d'^a (r, t) {3r^ - l) da (r, t) fooBore""* d 



+ 



r (r^ + 1) dr [r"^ + l)^ dr 



(r^ + l) 



— (ph {ar) + qKi (ar)) 



0. 



(53) 



An exact solution for a (r, t) for all r is not at our disposal, therefore, we shall search for solutions near r = and for 
r — oo. The solution for the homogenous part is 



<^H (r, t) = iCh {ar) + C2K, (ar)) . (54) 

Since a particular solution is not available we proceed to study the answers for a limited case, when r is small (to 
order r) 
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d'^a{r,t) d^a{r,t) lda{r,t) 



= 0, 



dt^ dr^ r dr 

a {r, t) = e-"*r {Cih {ar) + C2K1 (ar)) 



If we go to higher orders of r, (order for instance) we obtain 



dr 



Ar - 



da (r, t) 
dr 



boBor < -pa 



-Inl-a 



— hr r — 7 ) + pa 



whose solution can be expressed as 



2q 



= 



(55) 



(56) 



a{r,t) = e-'' 
in which 



CsWhitta kerM{--a^,-,2r^) + dWhitta ker M(- -a^ - , 2r'^' 
8 2 00 



1 2 1 
—a , — , 

8 ' 2' 



Y^r'^q (-8 + a^) In a - a^r'^q (-8 + a^) Inr + [9 (In 2 - 7) ~ p] a^r 
n = ( [(87 - 2 - 8 hi 2) - 2] g + ^r'^p) + IQq - 16r\. 

Here p, q, C3, C4, a and 7 are all constants and WhittakerM (— i ) stands for the Whittaker function^. For 

r — the function a (r, t) is finite. Now, to see the case when r goes to infinity we solve the differential equation 



Boboe- 



a{-8 + a^y 



4„2 



(S + n) (57) 



Of" 



Q^2 



0, (to order l/r) 



whose solution reads 



(58) 



(59) 



Here C5 and Cg are new integration constants. It can easily be checked that the ratio of this solution goes to zero for 
r — > 00 if we choose C5 = 0. 

IV. GEODESIC MOTION 

In this section we shall investigate the time-like (for d > 4) and null (for d = 3,4) geodesies by employing our line 
element given in Eq. (15). For d > 4 we divide the line element by dr (r is proper time) and for d = 3 by dA (A is an 
affine parameter) so that the Lagrangian can be expressed in the form 



L 



(d- 3)3^(r2 + 1)3^ 

\d\) \d\) 



\dT. 



,2 ( 



d-3 



( dr_Y _ V (^'' 
\ dT) \ dr 



(d-3)^(r 



5)^(r^ + l)^ \dT J 



d> 4 



d = 3 



(60) 



For the equations of motion with constant azimuthal angle (if ^constant for d > 4) and null geodesies for d = 3 we 
obtain 



dt 
dr 



dr 



fx = Hoe- 

dip Hi 

dX ~ 



do 

2 2 — 

2 

(d-3)<i-3 (r2 + l)3=3 



d> 4 



d = 3 



(61) 
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in which do , di^ Hq and Hi are all constants of integration. From the metric condition we find ^ (for d > A, ip 
constant) and (for d — 3) as follow 



dr I 

dr ~ ^ 



^ 2- ya^~id~3) "--^ (r2 + 1) <i-3 



al = dl-j:d^ 

1=1 



dX 

H2 = ^ 



d > 4 



d = 3 



(62) 



Hobo 



In effect, we obtain for d = A the relation between r and r but for d = 3 we want to find the relation between r and 



±(r-To) = (d-3)^ 
^ fco 



(r'^ + l)'^^ dr 
'fe2-(r2 + l) A 



(d-3)'i-3 



dr 



^0 J ry/r^-H^e"^'^ 

in which, tq and lySg £ire initial constants and we impose the restrictions so that 



d> 4 



d = 3 



(63) 
(64) 



< 1 



For d = 4 and 5, we have exact integrals given by 



T - To 



-^r^fc;^ - (r2 + 1)2 + V -('•0+ 



f7 



V2 



irVA:^-(r2 + l) + i (fc§ + l) arctan (^^|^;5^ 



in which 



n . _ 2.0) EUipticF (^^, ^-1 + 



(2A;o + 2) EUipticF 



d> 4 
d = 3 



d = 4 
d = 5 



1-1 



(fco + 1) 



(65) 



(66) 



Fig.s la and lb depict the behaviors of (63) (for d = 4) and (66) (for d = 5) , respectively. Now, we wish to consider 
the d = 4 null geodesies as well. The line element is 



ds^ - (r2 + 1)2 (dt2 - dr2 - dz') ""'^o 



(r2 + 1)- 



(67) 



Eq.s of motion imply for the affine parameter A 



dt 



an 



dz 



/3o 



dip _ (r^ + if 



(68) 



dA (r2 + l)2' dA (r2 + l)2' 

in which ao, and 7q are integration constants. We note that (67) and (68) correspond to Eq.s (1) and (6) of Q, 
respectively. From the null-metric condition ds2 = and upon shifting the independent variable to Lp we obtain 



(r2 + Ifdr 



(69) 
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where we have introduced the constant Aq as Aq = ■ Fig.s 2a and 2b are the plots corresponding to (64) and 

(69), respectively. 

Finally, we study the 4-dimensional time-like geodesies of a charged particle with charge e whose Lagrangian is 
given by 



The Euler-Lagrange equations of motion yield 



dip 

(r2 + l)2 Vd7 



eBobo 1 f dp 
2 (r2 + 1) [d^ 



(70) 



dt 



h 



dz 



Mo 



dp _l (-cTor^ - (To + hoai) (r^ + 1) 



(r2 + l)2' dT (r2 + l)2' dT 2 



by 



2 ' 



(71) 



for the integration constants Iq, Hq and ao- Since we shall be interested only in the r (r) behavior of the motion we 
derive the second order equation as follows 



-8r(r2 + 1) + 



4r 



(r2 + 1)3 



(r^ + 4r^ + 6r^ + Ar"^ + l) 



(72) 



+ 



2r36§(r2 + 1)3 



[ay° + Say + (2ag + 6^af ) - {2al + hi [-3a? + 2fj,l - 21^] ) 



r2 {3al-3albl)-al + albl 



0. 



For a set of chosen constants and initial values we plot the behavior of r (t) as depicted in Fig. 3. Our overall analysis 
shows that irrespective of the initial conditions r (r) ^ oo , with the increasing proper time. 



V. CONCLUSION 



We rederive the family of cylindrically symmetric magnetic universes in a particular metric ansatz which is con- 
formally flat on each constant azimuthal angle. These are non-black hole solutions where unlike their spherical 
counterparts the gravity of magnetic fields is not strong enough to make black holes. The energy conditions (in the 
Appendix) of the magnetic field are satisfied only in particular dimensions. Being inspired by the stability properties 
of the original 4-dimensional Melvin's magnetic universe and those of Gibbons and Wiltshire we prove also the sta- 
bility of present universes in a different gauge and in all dimensions including d = 3. Small radial perturbations of 
metric functions and the magnetic field (which automatically yields an electric field in accordance with the Maxwell 
equations) result in convergent expansions. Stability of the 3-dimcnsional case which was not considered in previous 
studies turns out to be weaker (i.e. convergence of perturbations are not satisfied simultaneously at r = and at 
r — )• oo). 

Geodesies show numerically that in running proper time uncharged particles are confined while null geodesies spiral 
around the center. Exact, particular geodesies in terms of the elementary functions are available in d = 5, whereas in 
d = 4 wo have elliptic functions. 

Acknowledgment: We wish to thank S. Habib Mazharimousavi for much valuable discussions. 
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Apendix A 
Energy conditions 

When a matter field couples to any system, energy conditions must be satisfied for physically acceptable solutions. 
We follow the steps as given infill ■ 
Weak Energy Condition (WEC): 



The WEC states that 



P > (Al) 

p + p, > 



In which p is the energy density and p^are the principle pressure given by 



p = n m 



The WEC conditions are trivially satisfied. 
Strong Energy Condition (SEC): 
This condition states that 



P+ > (A3) 

i=l 

P + Pi > Q 



For d=3, it means that 



> 0^3T°>0 (A4) 

1=1 

P + Pi > 0^ 2To° > 

which are satisfied. For 4 < d < 6 it is also satisfied because p + J^Pi = ^{'^ ^ 6)Tq . It can easily be seen also that 

i=l 

for 7 < d 

p+Ep. > o^-r°>o (A5) 

i=l 



P+P^ > 0^T°>0 



i.e. the SEC is violated. 



11 



Dominant Energy Condition (DEC): 

In accordance with DEC, the effective pressure should not be negative. This amounts to 



''» = ^%T:=^K>0 (A6) 



For having p^ff > o it is clear that d> 5. 
Causality Condition: 

In addition to the energy conditions one can impose the causality condition 



0<^<1 



(A7) 



This implies that 



< < 1 (A8) 



which is satisfied for d > 5. 



Figure captions 

Fig. 1 : Radial distance behavior as a function of proper time for timc-likc, neutral particle geodesies. From Eq. 
(63) in the text, the behavior of r (r) as the proper time r runs from zero to infinity. The starting points are chosen 
such that r = at T = 0, in both d = A (Fig. la) and d = 5 (Fig. lb)cascs. The fact that r is confined is clearly 
seen from these plots. This particular property is already implied from Eq. (65). Let us note that in both cases for 
simplicity we choose (+) sign and to — 0. Further, wc choose the constant fco = 2. 

Fig. 2 : Radial behavior as a function of the azimuthal angle for null geodesies in d = 3 (Fig 2a) and d = A 
dimensional (Fig. 2b) magnetic universes. The horizontal axis x {= r cos ip) and vertical axis y {= r sin <p) are plotted 
numerically in each case from the expressions given in Eq. (64) and Eq. (69) . For simplicity wc choose the constants 
H2 = 0.01, 60 = 0.01 , Hq = 1, c = I and Aq = 1. Let us note that we have chosen the (+) sign in both (64) and (69) 
, which give outward orbits around the center. Obviously, the choice (— ) should yield inward orbits. 

Fig. 3 : This is a numerical plot of the intricate second order differential equation (72) in d = 4. For technical 
reason we choose the constants ai, bo, lo, A*o ^^'^ ^0 all equal to one. Two different initial conditions are displayed 
(A and B) which reveal the pattern of increasing radial distance in proper time. 
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